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Purpose: This study aimed to investigate useful parameters for estimating gastroc- 
nemius (GCM) muscle volume (MV) using ultrasonography (US) and anthropom- 
etry in children with spastic cerebral palsy (CP). Materials and Methods: Eigh- 
teen legs from nine children with spastic CP aged 2 to 6 years were investigated in 
this study. Tibial length (TL) of each leg was measured and muscle thickness 
(MT) and anatomical cross-sectional area (aCSA) of GCM muscles were assessed 
using US. The volume of the GCM was measured by magnetic resonance imaging 
(MRI) scans. The relationship of TL, MT, and aCSA with MV measured by MRI 
was investigated. Simple and multiple regression analyses were performed to es- 
tablish muscle volume prediction equations. Results: Resting MT, aCSA, and TL 
were highly related to MV of both medial and lateral head of GCM determined by 
MRI. The MV prediction equation based on simple regression analysis resulted in 
r 2 values ranging from 0.591 to 0.832 (p<0.05). The r 2 values were higher using 
aCSA as independent variable than using MT. The MV prediction equation based 
on multiple regression analysis resulted in r 2 values ranging from 0.779 to 0.903 
(p<0.05). However, the relatively high standard error of the estimate values ranged 
from 18.0-33.6% on simple regression and 15.5-25.6% on multiple regression. 
The contribution of aCSA was higher than that of MT for predicting MV of GCM. 
Conclusion: Our study demonstrated the suitability of US assessment of aCSA 
and MT combined with TL for estimating MV of GCM in children with spastic 
CP and showed that aCSA is more useful parameter than MT. 
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Cerebral palsy (CP) describes a group of disorders affecting the development of 
movement and posture, and thus limiting activity, that are attributed to non-pro- 
gressive disturbances in the developing fetal or infant brain. 1 According to pub- 
lished literature, there is consistent evidence for a reduction of muscle volume 
(MV) in paretic limbs of children with CP, compared with both non-paretic limbs 
of children with CP and limbs of typically developing peers. 2 Deficits in MV in 
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children with spastic CP may increase as they grow older. 3 

MV is a good predictor of maximal force generating ca- 
pacity, 49 and accordingly, the lack of growth in muscle vol- 
ume may contribute to muscle weakness, which may in 
turn be related to functional limitations in these children. 1011 
A previous report suggested that measurement of muscle 
size may be an alternative measure of quantitative muscle 
evaluation for people with severe CP in whom direct mea- 
surement of muscle strength is difficult. 11 In this context, 
measurement of MV in a clinical setting would seem to be 
useful as a means of selection for intervention and for de- 
lineating the consequences of interventions, such as botuli- 
num toxin injection or strengthening exercises, on muscle 
volume. 1213 

Previous attempts were made to estimate MV using an- 
thropometric data, 1415 as well as by magnetic resonance im- 
aging (MRI), computed tomography (CT), and ultrasonogra- 
phy (US); CT and MRI are considered the "gold standards" 
of MV measurement. 1617 More recently, three-dimensional 
(3D) US was shown to be a valid and reliable measurement 
of muscle volume. 18 20 However, two-dimensional bright- 
ness mode (2D B-mode) US can be used to produce high- 
quality images of muscle morphology, and has the same ad- 
vantages as MRI in distinguishing between muscle and fat 
tissues. 17 Moreover, 2D B-mode US is simple and easy to 
perform and requires less time for analysis, compared with 
MRI, CT, and 3D US. Consequently, 2D B-mode US is 
commonly used to measure muscle size, such as anatomical 
cross-sectional area (aCSA) and muscle thickness (MT). In 
the published literature, measurement of aCSA and MT us- 
ing 2D US exhibited excellent reliability in evaluating chil- 
dren and adolescents with CP. 2122 Additionally, there have 
been some attempts to estimate MV based on parameters 
obtained from 2D US scan images. 23 27 MV prediction equa- 
tions have been developed based on MT measurement us- 
ing US and measured tibial length against muscle weight of 
cadaveric human legs 23,24 or MV measured by MRI in 
healthy adults. 25 27 The results of such studies suggest that 
MT measurement using US is reliable and useful for esti- 
mating MV 25 " 28 

Anatomical CSA of a muscle is one of several muscle ar- 
chitecture parameters that are closely related to maximum 
muscle strength. According to Blazevich, et al., 6 aCSA is 
the second-best predictor of maximal muscle torque. How- 
ever, MV estimation using aCSA measured from a single 
US scan image has not yet been investigated in children 
with cerebral palsy. 



Equinus foot is the most common deformity in children 
with spastic CP, and thus, the gastrocnemius (GCM) mus- 
cle is the most common target for botulinum toxin injec- 
tion. The optimal dose of the toxin may be related to MV. 
Recently, muscle architectural changes after toxin injection 
have become an interesting issue with respect to the influ- 
ence of the toxin on growth in muscle volume. 29 31 Never- 
theless, the suitability of 2D B-mode US for MV estimation 
of the GCM in children with spastic CP has not been inves- 
tigated. Therefore, the aim of this study was to investigate 
the application of US and identify useful parameters for the 
prediction of GCM MV in children with spastic CP. 



MATERIALS AND METHODS 



Participants 

Children with spastic CP who were admitted or visited our 
out-patient clinic from January 2010 to May 2012 were re- 
cruited for this study upon fulfilling the following inclusion 
criteria: 1) children with spastic cerebral palsy, 2) Gross 
Motor Function Classification System (GMFCS) level I to 
IV, and 3) age of 2 to 6 years. The exclusion criteria were as 
follows: 1) previous history of GCM muscle trauma, 2) 
combined with muscular disease, 3) previous history of sur- 
gical intervention of lower extremity, 4) any metal implant 
in the body, or 5) chemodenervation therapy within 6 
months. 

Ethical approval was granted by the Institutional Review 
Board and ethics committee of Severance Hospital (# 4-2009- 
0650). As all the children in this study were younger than 
18 years, informed consent was obtained from their parents 
for participation in the study. 

Ultrasound measurements 

US images of the GCM were taken by a physiatrist using 
2D B-mode and real-time ultrasonography (ACCUVLX VI 0 
system, Samsung Medison Co, Seoul, Korea) with a linear- 
array probe (5 to 12 MHz). 

The subjects were positioned prone on the examination 
table with their feet hanging over the edge of the table. US 
images over the muscle belly of both medial and lateral 
heads of GCM muscles were taken in a resting ankle posi- 
tion. For image acquisition, the ultrasound probe was trans- 
versely positioned on the skin with minimal compression. 
US images were obtained following the recommendations 
of Benard, et al. 32 in order to minimise measurement error. 
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The ultrasound probe was centred on the mid-longitudinal 
axis, at 25% of the tibial length from the popliteal fossa. Ul- 
trasound images were taken over a range of tilt and rotation 
angles of the ultrasound probe at 5° increments. Tilt and ro- 
tation are expressed as deviations from the longitudinal ori- 
entation perpendicular to the skin for all conditions. This 
was performed three times for each muscle by the same in- 
vestigator to obtain images of maximal MT. 

MT was defined as the longest distance between the up- 
per muscular fascia and the lower muscular fascia visible 
on the image. Anatomical CSA represents the cross-section 
perpendicular to the longitudinal axis of a muscle, and is 
defined as the area surrounded by the upper muscular fas- 
cia, the lower muscular fascia, and the intramuscular sep- 
tum. 33 We calculated anatomical CSA using the same image 
of maximal MT. Fig. 1 shows an example of imaging mea- 
surement of MT and aCSA. 

GCM volume measurement by MRI 

The entire GCM muscle was imaged using axial, gapless, 
Tl -weighted sequences. MR images were viewed and pro- 
cessed using OsiriX® software (version 3.6.1, University 
Hospital of Geneva, Switzerland, 32-bit, http://www.osirix- 
viewer.com). The MRI images were obtained in the resting 
ankle position. Volumetry of the GCM muscle was per- 
formed manually. A physiatrist manually outlined the junc- 
tion of the GCM muscle fascia and intermuscular septum, 
first cranially and then caudally, in the axial image stack. 
OsiriX® was used to calculate the GCM muscle volume by 
multiplying all outlined areas of the slices by the slice thick- 
ness. 8 

Statistical analysis 

Univariate analysis was performed to determine the corre- 
lation between muscle volume of medial and lateral heads 
of the GCM and the TL, MT, or aCS A using Pearson corre- 
lation analysis. Simple and multiple regression analyses were 
performed using muscle volume measured by MRI as the 
dependent variable and [ttx(MT/2) 2 xTL] or (CSAxTL) as 
the independent variables. Multiple regression analyses 
were also performed using TL, MT, and aCSA as indepen- 
dent variables. Using the established muscle-volume pre- 
diction equations, the muscle volume for medial and lateral 
heads of the GCM muscle was calculated for each leg. To 
assess the criterion validity of the equations in predicting 
the measured muscle volume, the standard error of the esti- 
mate (SEE) and validity correlation (r 2 ) were calculated us- 



ing the measured muscle volume as the criterion. In multi- 
ple regression analysis, each independent variable was 
calculated as an index representing its relative contribution 
to estimating the measured muscle volume. Contribution 
was calculated by partial r 2 . The SEE was calculated using 
the equation SD diff Hi, where SD diff is the standard de- 
viation of the difference scores between the estimated and 
measured muscle volume. The SEEs were expressed as ab- 
solute values and values relative to the measured muscle 
volume. 



RESULTS 



Characteristics of participants 

Nine children with bilateral involvement spastic CP (six boys 
and three girls; mean age±SD, 3.8±1.5 years, range, 2-6 
years) were included in this study. The distribution of GM- 
FCS was as follows: Level I (n=5), Level II (n=l), Level III 
(n=l), and Level IV (n=2). Only one patient had received 
BoNT-A injection within 6 months before the study. 

In total, 18 lower limbs were evaluated. For anthropomet- 
ric measurements, TL was measured from the lateral knee- 
joint line to the most prominent point of the lateral malleo- 
lus. Mean TL was 19.9±2.7 cm (range 17.1-25.0 cm). 

Muscle volume, thickness and cross-sectional area 

The MV of GCM measured by MRI was 19.70±9.29 cm 3 
(range; 7.18-38.95 cm 3 ) for the medial head and 11.92± 
9.12 cm 3 (range; 4.14-28.05 cm 3 ) for the lateral head. The 
MT of GCM muscles measured using US was 0.99±0.27 
cm (range; 0.48-1.31 cm) for medial head and 0.82±0.20 




Fig. 1. Ultrasonographic measurement of muscle thickness and anatomical 
cross-sectional area. Muscle thickness (MT): the longest distance between 
the upper muscular fascia and the lower muscular fascia. Anatomical cross- 
sectional area (aCSA): the area surrounded by the upper muscular fascia, 
the lower muscular fascia, and the intramuscular septum. 
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cm (range; 0.46-1.17 cm) for lateral head. The aCSA of 
GCM muscles was 3.05±1.03 cm 2 (range; 1.35-4.72 cm 2 ) 
for medial head and 1.92±0.73 cm 2 (range; 0.76-2.98 cm 2 ) 
for lateral head (Table 1). 



SEE values of the equations were 3.5 cm 2 (18.0%) and 2.6 
cm 2 (21.6%) (Table 2). The relationship between estimated 
and measured MV using simple regression analysis of chil- 
dren with cerebral palsy is shown in Fig. 2. 



Univariate analysis 

Tibial length, MT and aCSA from US imaging, were highly 
related to MV determined by MRI for both the medial 
heads (r=0.793, 0.659, 0.810, respectively, /?<0.05) and lat- 
eral heads (r=0.769, 0.621, 0.824, respectively, /K0.05) of 
GCM muscles. 

Simple regression analysis 

In simple regression analyses, [7ix(MT/2) 2 xTL] and (CSAx 
TL) were significantly correlated with the measured MV 
for both the medial and lateral heads. Simple regression 
equations using MT as an independent variable, based on 
the relationship between MT and measured MV, resulted in 
values of r 2 =0.653 (p<0.05) for medial heads and r 2 =0.591 
(p<0.05) for lateral heads. The SEE values in the equations 
were 5.6 cm 2 (28.6%) for medial heads and 4.0 cm 2 (33.8%) 
for lateral heads. The MV prediction equation using aCSA 
as independent variable resulted in r 2 =0.862 (p<0.05) for 
medial heads and r 2 =0.832 (p<0.05) for lateral heads. The 

Table 1. Characteristics of Patients 



Mean±SD (range) 


Sex (M:F) 


6:3 


Age (months) 


49.7±17.1 (24-74) 


Weight (kg) 


15.7±4.9 (10.9-28.0) 


Tibial length (cm) 


19.9±2.7 (17.1-25.0) 


GCM muscle volume (cm 3 ) 


Medial 


19.70±9.29 (7.18-38.95) 


Lateral 


11.92±9. 12 (4.14-28.05) 


Muscle thickness (cm) 


Medial 


0.99±0.27 (0.48-1.31) 


Lateral 


0.82±0.20 (0.46-1.17) 


Anatomical cross-sectional area (cm 2 ) 


Medial 


3.05±1. 03 (1.35-4.72) 


Lateral 


1.92±0.73 (0.76-2.98) 



Multiple regression analysis 

Multiple regression analysis revealed that MT, aCSA, and 
TL were significant contributors to the prediction of MV 
measured by MRI. The muscle volume prediction equa- 
tions based on multiple regression analyses are given in Ta- 
ble 3. The final equation based on multiple regression analy- 
sis using two independent variables of TL and MT resulted in 
r 2 =0.831 (/?<0.05) for medial heads and r 2 =0.779 (p<0.05) 
for lateral heads. Based on this equation, the relative contri- 
bution of MT to the variation in MV measured by MRI was 
20.1% for the medial head and 18.8% for the lateral head. 

The final equation based on multiple regression analysis 
using two independent variables of TL and aCSA resulted in 
r 2 =0.903 (/?<0.05) for medial heads and r 2 =0.858 (p<0.05) 
for lateral heads. Based on this equation, the relative contri- 
bution of aCSA to the variation in MV measured by MRI 
was 65.6% for the medial head and 67.8% for the lateral 
head (Table 3). 



DISCUSSION 



In published literature, MV in children with CP was mea- 
sured using MRI 34 36 and 3-D US. 312 ' 2137 Using these tech- 
niques MV can be reliably and validly measured in children 
with CR On the other hand, 2D US allows for excellent dif- 
ferentiation between the connective tissue, fat, and muscle 
tissue, and is also easier to perform and less time-consuming 
for analysis. Therefore, 2D US is commonly used as a means 
of assessing muscle size. Previous studies showed that MT 
measurements obtained using 2D US are useful for estimat- 
ing MV in different muscle groups. 24,26,27 To our knowledge, 
two studies have tried to estimate MV of ankle plantar flex- 
or from a single 2D US image. 24,27 One study investigated 
the suitability of MT measurement by US for estimation of 



GCM, gastrocnemius. 

Table 2. Simple Regression Equations for Predicting Muscle Volume in Children with Cerebral Palsy 





Equation 


r 2 (p value) 


SEE (cm 3 , %) 


GCM medial head 


y=0.868xTLx7i(MT/2) 2 +5. 112 


0.653 (0.000) 


5.6 (28.6) 




y=0.332xTLx a CSA -0.958 


0.862 (0.000) 


3.5(18.0) 


GCM lateral head 


y=0.802xTLx7i(MT/2) 2 +2.641 


0.591 (0.000) 


4.0 (33.8) 




y=0.298xTLx a CSA+0.203 


0.832 (0.000) 


2.6 (21.6) 



TL, tibial length; MT, muscle thickness; aCSA, anatomical cross-sectional area; SEE, standard error of the estimate; GCM, gastrocnemius. 
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Fig. 2. Relationship between estimated (based on ultrasound and arthropometry) and measured (MRI) muscle volume for gas- 
trocnemius muscle, using simple and multiple regression analyses. (A) y^O.^xJlxniMJUf+bAM, r 2 =0.653 (p=0.000), SEE= 
5.645 cm 3 (28.6%). (B) yp0.332xTLxCSA-0.958, r 2 =0.862 (p=0.000), SEE=3.559 cm 3 (18.0%). (C) y 2 =0.802xTLxn(MT/2) 2 +2.641, r 2 =0.591 
(0.000), SEE=4.1 cm 3 (33.8%). (D) y 2 =0.298xTLxCSA+0.203, r 2 =0.832 (p=0.000), SEE=2.6 (21.6%). (E) Y 1= 2.271 xTL+!5.982xMT-41 .493, 
r 2 =0.831 (0.001), SEE=4.1 (20.6%). (F) Y 1= l .896xTL+5.192xCSA-34.050, r 2 =0.903 (0.000), SEE=3.1 (15.5%). (G) Y 2 = 
1.479xTL+13.347xMT-28.676, r 2 =0.779 (0.003), SEE=3.1 (25.6%). (H) Y 2 =1.101xTL+4.903xCSA-19.516, r 2 =0.858 (0.001), SEE=2.4 
(20.5%). TL, tibial length; MT, muscle thickness; CSA, cross-sectional area; SEE, standard error of the estimate. 
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Table 3. Multiple Regression Model for Predicting Muscle Volume in Children with Cerebral Palsy 



Muscle group 


Variable 


r 2 (p value) 


SEE (cm 3 , %) 


Regression 
coefficient 


Standard regression 
coefficient 


Contribution % 


CjCM medial 


TL 
MT 


0.831 (0.001) 


4.1 (20.6) 


2.271 

1 J.70Z, 


0.658 


62.9 


head 


TL 

aCSA 


0.903 (0.000) 


3.1(15.5) 


1.896 
5.192 


0.549 
0.578 


24.8 
65.6 


GCM lateral 


TL 
MT 


0.779 (0.003) 


3.1 (25.6) 


1.479 
13.347 


0.651 
0.449 


59.1 
18.8 


head 


TL 

aCSA 


0.858 (0.001) 


2.4 (20.5) 


1.101 
4.903 


0.484 
0.590 


18.0 
67.8 



TL, tibial length; MT, muscle thickness; aCSA, anatomical cross-sectional area; SEE, standard error of the estimate; GCM, gastrocnemius. 



MV of ankle plantar flexor against MV measured by water 
displacement of dissected muscles in cadaveric human legs, 
and the regression between estimated MV and measured MV 
was r 2 =0.373 by a simple regression equation and r 2 =0.497 
by a multiple regression equation. However, the results ob- 
tained in cadavers may differ from data obtained from liv- 
ing humans because of cell shrinkage resulting from freez- 
ing of human tissue, which might account for the discrepancy 
in their data. On the other hand, Miyatani, et al. 27 investigat- 
ed the accuracy of estimating ankle plantar flexor MV in 
healthy adults using MT measurement from a single US 
image, and demonstrated that the regression between MV 
estimated by a MV prediction equation and that measured 
by MRI was r 2 =0.649 by a simple regression equation and 
r 2 =0.820 by a multiple regression equation. 27 They conclud- 
ed that MT on US images is useful for estimating ankle 
plantar flexor MV The results of our study are in accor- 
dance with those of Miyatani, et al., 27 and confirmed that 
MT on 2D US is a useful parameter when combined with 
TL for estimating GCM MV in children with spastic CP. 

The functional importance of changes in the aCSA of a 
muscle can be seen in its contribution to total force produc- 
tion, 38 and aCSA is therefore often analysed in studies of 
strength training, in which it shows a strong relationship 
with the maximum capacity of muscle force production. 4 ' 39 40 
In the literature, aCSA measured using 2D B-mode US scan- 
ning shows a high correlation with other imaging techniques 
such as MRI or CT. 16 25 41 Therefore, US has been common- 
ly used to measure human skeletal muscle aCSA in vivo in 
different subjects and different muscles. 25 ' 4042 For children 
and adolescents with CP, aCSA measurement using 2D US 
had excellent reliability 21 Using this technique, aCSA of large 
muscles can be reliably and validly measured using 2D US. 17 
We measured aCSA of each head of GCM muscles from 
routine single transverse plane US images without use of the 



panoramic scanning technique because the aCSA of each 
head of the GCM muscles can be shown in whole in young 
children. In multiple regression equations using measured 
aCSA and TL, the r 2 values were 0.903 for medial heads 
and 0.858 for lateral heads. These high values of r 2 suggest 
the usefulness of aCSA measurement for predicting MV. 
According to a previous study, estimations of MV from the 
aCSA of the tibialis anterior measured by 1 1 axial US scans 
were highly reliable and valid when compared with MV 
measured by MRI. 25 In addition, the study of Morse, et al. 43 
demonstrated that aCSA from a single MRI scan was highly 
correlated with the MV of quadriceps measured by 1 1 con- 
tinuous transverse plane MRI scans along the entire length of 
the femur. These findings suggest that aCSA may be useful 
for estimating MV Compared with MRI and the multiple 
axial US scans used in previous studies, aCSA measure- 
ments from a single US image are easier to perform and ac- 
cess, and also require much less time for analysis. In this 
context, the clinical implication of the MV prediction equa- 
tion obtained in our study is noteworthy. As far as we are 
aware, this is the first study to show the suitability of aCSA 
measurements using 2D US for estimating MV in children 
with CP. 

Nevertheless, the relatively high SEE values, ranging from 
18.0-33.6% in simple regression and 15.5-25.6% in multi- 
ple regression, suggest a limitation in the MV prediction 
equation for use in individual assessment in longitudinal 
studies. To reduce the SEE values, increasing the number of 
levels of measurement and standardisation of US imaging 
with respect to the ankle joint position, the severity of spas- 
ticity and ankle joint limitations in motion should be taken 
into consideration. Further studies are needed to find the 
most efficient and accurate way to estimate MV for individ- 
ual assessment in longitudinal studies. 

In conclusion, MV prediction equations using MT and 
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aCSA of GCM muscles measured by 2D US combined with 
TL were obtained based on MV of GCM measured by MRI 
in children with spastic CP. The high values of r 2 in each 
multiple regression MV prediction equation obtained in our 
study indicate that MV estimates using MT or aCSA mea- 
surement using 2D US are suitable when combined with 
tibial length. In addition, aCSA is more useful than MT for 
estimating the MV of GCM. The final MV prediction equa- 
tions obtained in our study may provide a quick and easy 
way to estimate the volume of GCM muscles in these chil- 
dren, which may ultimately be useful for population studies 
or field surveys. However, at present the SEE values ob- 
served in our study indicate that it might be better to limit 
use of these prediction equations to individual assessment 
in longitudinal studies. Further attempts to improve the ac- 
curacy of MV estimation using 2D US are needed. 
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